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Abstract 

We introduce a Wilson line distribution function W T (v) to study gluon sat- 
uration at small Feynman x F , or large r = ln(l/x F ). This new distribution 
can be obtained from the distribution W T (a) of the Color Glass Condensate 
model and the JIMWLK renormalization group equation. W T (v) is physically 
more relevant, and mathematically simpler to deal with because of unitarity 
of the Wilson line v. A JIMWLK equation is derived for W T (v); its proper- 
ties are studied. These properties are used to complete Mueller's derivation of 
the JIMWLK equation, though for W T (v) and not W T (a). They are used to 
derive a generalized Balitsky-Kovchegov equation for higher multipole ampli- 
tudes. They are also used to compute the unintegrated gluon distribution at 
x F = 0, yielding a completely flat spectrum in transverse momentum squared 
k 2 , with a known height. This is similar but not identical to the mean field 
result at small k 2 . 



I. INTRODUCTION 

Soft gluons are produced by radiation from more energetic partons. Since the number of 
sources increases at small Feynman x F , the soft gluon density x F G per unit rapidity interval 
increases with r = ln(l/x F ). In fact, both the DGLAP equation |J and the BFKL equation 
0] predict a growth so fast that the unitarity bound r 2 is violated. To restore unitarity, a 
new mechanism is required to slow down the growth towards x F = 0. The momentum 
Q s at which this mechanism sets in is known as the saturation momentum. 

The phenomenological implications for the presence of a saturation momentum have 
been discussed in many papers [|J], but it is not yet clear whether saturation has been 
observed experimentally. On the theoretical side, the important thing to note for our present 
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discussion is that soft gluons can be treated as a classical color potential a(x), because of 
its large density at small x F . In this background, energetic partons interact with soft gluons 
through their Wilson lines. It is this interaction that is responsible for the saturation process. 

With the replacement of soft gluons by a classical Yang-Mills potential, the growth of soft 
gluon density is determined by the r dependence of the distribution function W T (a). This 
dependence is given by a renormalization group equation known as the JIMWLK equation 

In this paper we introduce and study a distribution W T (v) for the Wilson line v. We 
will show that this distribution can be obtained from W T (a), and that it still satisfies the 
JIMWLK equation. Compared to W T (a), it has the advantage of being more directly phys- 
ically relevant, because many physical processes can be described in terms of Wilson lines 
or dipole amplitudes. Moreover, the Wilson lines v are unitary matrices, living on the color 
group manifold which is compact. This compactness brings with it a number of mathe- 
matical advantages not shared by W T (a), whose argument a lives in a non-compact linear 
space. 

W T (v) can be used to derive a number of properties not easily obtainable directly from 
W T (a). We will use these properties to complete Mueller's proof of the JIMWLK equation 
H, though the proof is valid for W T (v) and not for W T (a). They will be used to derive 
a generalized Balitsky-Kovchegov (BK) equation for multipole amplitudes. We find that 
once the non-linear BK equation for the dipole amplitude is solved, all subsequent multipole 
amplitudes can be obtained by solving only linear equations. This means saturation of the 
dipole amplitude automatically leads to saturation of higher multipole amplitudes. 

W T (v ) will also be used to compute the asymptotic behavior of unintegrated gluon dis- 
tribution at x F =0. We get a flat distribution in the gluon transverse momentum k 2 , with 
a height given by eq. ( |6.7j ). This is to be contrasted with the mean field result ( |2.1(j| ) which 
yields a logarithmic dependence on k 2 with an undetermined normalization. Pre-asymptotic 
corrections will also be breifly discussed. 

In the next section we start with a short review of saturation, the JIMWLK equation, 
the related BK equation, the BFKL equation, as well as some of their solutions. In Sec. 3, 
the Wilson line distribution W T (v ) is introduced, together with some mathematical prelim- 
inaries. The JIMWLK equation for W T (v ) is derived, and its properties studied. In Sec. 4, 
the missing steps of Mueller's derivation of the JIMWLK equation are supplied. In Sec. 5, 
the infrared divergence encountered in the JIMWLK equation is discussed. It is shown that 
certain multipole functions which we call multipole traces are free of these divergences. A 
generalized BK equation is derived to describe the evolution of the multipole amplitudes. 
In Sec. 6, we discuss the property and implications of the asymptotic solution of W T (v) at 
x F = 0. In particular, the unintegrated momentum spectrum of the gluon density is derived. 
Correction to the asymptotic limit when x F ^ is briefly discussed. Certain mathematical 
details can be found in Appendices A and B at the end. 

After this paper was submitted for publication, we were informed by the Referee of an 
interesting recent paper 0] in which W T (v ) was also introduced and its JIMWLK equation 
derived. It went on to give a random-walk interpretation of the JIMWLK equations, with 
r playing the role of time, v the position in the group manifold, and a the velocity. 
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II. A BRIEF REVIEW 



The number of gluons emitted by a valence quark, per unit rapidity r = \n(l/x F ) and 
per unit transverse momentum squared k 2 , is given in perturbation theory by a s CV/7rk 2 , 
where a s = g 2 /47r is the QCD fine structure constant. Cp is the Casimir number in the 
fundamental representation, being (N 2 — 1)/2N C for the gauge group SU(N C ), and \N C for 
U(N C ). A nucleus with atomic number A has AN C valence quarks, so its unintegrated gluon 
distribution is dN/drdk 2 = d(x F G)/dk 2 = a s C F (AN c ) /irk 2 . 

When integrated, this formula encounters an infrared divergence at small k 2 , brought 
about by the long range gluon field of the unshielded valence quarks. However, quarks are 
confined inside color-singlet nucleons, so such a long range force is absent beyond the nucleon 
radius a. Thus x F G(x F ,Q 2 ) = (a s Cp(AN c )/Tr)ln(Q 2 a 2 ). 

In the central rapidity region where Feynman x F is small, the gluon density is much 
larger than the amount given by the perturbation formula. This is so because soft gluons 
can be radiated also from energetic gluons and sea quarks, not just the valence quarks 
considered so far. According to the DGLAP equation, the soft gluon density grows like 
exp(/« v /r) for some positive constant k, and, according to the BFKL equation 0, it grows 
like exp(4a s ln(2)A c r/7r). Both exceed the unitarity limit r 2 , so a new mechanism must 
kick in to dampen the growth and restore unitarity at small x F [[|. This effect is known as 
saturation. 

Saturation is thought to arise from a non-linear mechanism which occurs when gluons 
are sufficiently dense to interact among themselves []8]||. The number of gluons per unit 
rapidity interval is x F G. In a nucleus of radius Ra, the transverse area per gluon is therefore 
nR 2 A /x F G(x F) Q 2 ). The average color-charge squared of a gluon is N C /(N 2 — 1), their inter- 
action strength is ~ <y s /ir, and their natural size is ~ Hence the cross section for two 
gluons to interact can be estimated to be (a s /vr)(A c /(A 2 - 1))(tt/Q 2 ). If the cross section 
is larger than the transverse area per gluon, then interaction will take place to set off the 
non-linear mechanism. The onset therefore occurs at a momentum scale Q s such that 

rrR 2 A /x F G(x F ,Q 2 ) = c(a s /n)(N c /(N 2 - l))(7r/Q 2 ), (2.1) 

or equivalently, 

a 

A 



' ** ■ <2 ' 2) 



A constant c has been inserted to account for the qualitative nature of this argument. Even 
when Q 2 is obtained from a detailed calculation, the constant c is still somewhat ambiguous 
because transition into saturation does not occur sharply. Thus one finds a number of c's 
used in the literature. For example, c = 1 in ||10| , c = tt 2 in when estimated from the 



mean field approximation on the large k side, and c = IQti C\ when estimated on the small 
k 2 side, where C\ is some unknown constant. And, c = ir in [TTI . 



Q s can also be defined through the unintegrated gluon density d(x G G)/dk 2 . When we 
reduce k 2 from infinity, this density increases until a point k 2 = Q 2 when gluons become 
sufficiently dense to set off the non-linear mechanism. From there on we enter a saturation 
region with a much slower growth. However, this definition is also ambiguous unless the 
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slowdown occurs fairly sharply, which turns out to be the case at x F = 0. As we shall see in 
Sec. 6A, at x F = 0, the saturation region is large and the unintegrated spectrum d(x F G)/dk 2 
in this region is absolutely flat in k 2 . This then allows c to be determined unambiguously 
to be c = 87r 3 . 

Using the BFKL solution for x F G as a qualitative estimate, and assuming that x F G is 
proportional to A, we see from (2.2) that Q 2 grows with a power of l/x F and A 1 / 3 , making 



F 

it large for large nucleus or small x F . ( |2.2[ ) also implies that the gluon number per unit 
transverse area at saturation is ~ Q 2 (x)/a s . 

The large number of gluons present at saturation allows them to be treated as a classical 
Yang-Mills (YM) potential a a (x). The superscript a is the color index, and x = (x~,x) are 
the lightcone (LC) coordinates, defined for a hadron moving along the +z direction to be 
x ± = (t ± z)/y/2 and x = (x 1 ,^ 2 ). It is also convenient to introduce the spacetime rapidity 
variable y = ln(a; _ P + ), where P + is the + component of the hadron momentum, and the 
gluon potential a^(x) = x~a(x~ ,x). 

For a fast moving hadron (or nucleus), Lorentz contraction forces a(x) to be concentrated 
around x~ = 0, and time dilation makes it effectively (LC) time {x + ) independent. The soft 
gluons are produced by partons within the hadron, so one can assume o^(x) = for y > r 

!• 

Energetic partons, whether in the same hadron or not, interact with the soft gluons 
through the Wilson-line factors [] 



i>t(x) = Pexp y+ig J dx a a (x)t a ^j = Pexp (+ig J dy o^(x)t a 
v (x) = Pexp (—ig J dx~a a (x)ta^j = Pexp ^— ig J dy a^(x)t a ^ , 
V f (x) = Pexp f^+ig J°° dx-a a (x)T^j = Pexp (+ig J T dy a a y (^)T^j , 
V(x) = Pexp (-ig J°°^ dx-a a (x)T a ^j = Pexp (-ig £ dy a*(x)T a ) . (2.3) 

where P and P indicate respectively path-ordering and anti path-ordering. The first two 
expressions describe the propagation of quarks and anti-quarks, respectively, through the 
dense background of the soft gluons, and the last two expressions describe the propagation 
of gluons. These Wilson lines play a central role in the rest of the paper. 



We shall use upper case letters to denote the adjoint representation and lower case letters to 
denote the defining representation. In this notation, the generators in the defining representation 
will be denoted by t a , and they will be normalized to be tr(£ a t&) = \5 a b- The generators in the 
adjoint representation are denoted by T a . They are related to the totally antisymmetric structure 
constants by (T a )b c = ifbac- Hence (T a )i, c is imaginary and totally antisymmetric in the three 
indices. Similarly, the quark and anti-quark Wilson lines will be denoted by the lower-case letters 

and v, and the gluon Wilson line will be denoted by the upper-case letters and V. A 
slight drawback of this convention is that we are forced to denote the dipole amplitude ( |2.8p by a 
lower-case letter s T , whereas the usual notation for it is S T . 
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In this representation of soft gluons by a classical background field, gluon distribution is 
determined by the distribution W T (a) of the YM potential. W T depends on r because the 
number of sources available to emit soft gluons increases at small x F . 

The resulting change of the distribution functional W T (a) can be shown to satisfy the 
JIMWLK renormalization group equation || 



dW T (a) 
dr 



-HW T (a), 



(2.4) 



where 



H= - I rf 2 zrf 2 xrf 2 y^(xy|z)C(xy|z) 

7T J 



0(xy|z) 



s 



0(x)-W(z)l [V(y)-V(z)} 



S 



cb 



K(y)' 



and 



K{xy\z) 



1 (x 



47r 3 (x - z) 2 (y - z) 2 



(2.5) 



(2.6) 



A consequence of Q2.4Q and (|2.5| ) is that the normalization / X>[a]W T (a) is independent 
of r. We will normalize it to be 1, so that the average of any functional of a is given simply 
by (F) T = fV[a]F(a)W T (a). 

The functional derivatives 5/5a^ of V' and V in ( |2.5[ ), and similarly of v* and v that we 
will encounter later, are 

gM= WV t ( . ) , (x _. ) , 



fat(z) 
5a°(x) 

5a? (x) 
<5i>(z) 
5a? (x) 



iflrt a v t (z)£(x - z), 
-i5fV(z)T a 5(x - z), 
-^V^(z)t a (5(x - z). 



(2.7) 



For calculational simplicity it is useful to note that V^(x) — V\z) in fl2.5|) may be 
put in front of the operator S/Sa^(x). This follows from (|2.7| ) and the observation that 

(Fa) ac ifaac 0. 

A particularly important physical quantity to study is the dipole amplitude fi~2| 



Sr(x,y) = — (tr [u t (x)i;(y)]) T . 



(2.8) 



At coincident points 



s T (x,x) = 1 



(2.9) 
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because of unitarity of v. It can be shown form ( |2.4j) that the dipole amplitude satisfies the 
Balitsky equation fl3 



<9s T (x,y) 
dr 



a s N c 
' 2tt 2 



d 2 z- 



lx-yj 



x — z) 2 (y — z) 2 



s r(x,y) - -^(tr 



vHx)v(z) 



tr 



(2.10) 



Note that the infrared divergence occuring at large z in fl2.5|) is absent in the Balitsky 
equation. 



For large N c , the last term factorizes and we arrive at the Kovchegov equation ||14 
<9s r (x,y) _ a * N c [ d 2 z (x-y) 2 



dr 



2vr 2 



x — z) 2 (y — z) 2 



Sr(x,y) - s T (x,z)s r (z,y) \. 



(2-11) 



Since s T (x, x) = 1, the quantity t T (x, y) = 1 — s r (x, y) is expected to be small when 
x ~ y. If we deal with soft gluons outside of the saturation region, when a(x) is small, the 
Wilson lines v and are close to 1 anyway, so we expect to be able to drop the quadratic 
term t T (x, z)t T (z, y). In this way we get the dipole form of the BFKL equation [Q] 



dtr(x,y) 
dr 



a s N c 
' 2tt 2 



d 2 z- 



z) 2 (y 



t T (x, y) - t T (x, z) - t T (z, y) 



(2-12) 



For |x — y| 2 ^> 1/Q 2 and inside the saturation region, the strong classical YM potential 
a causes large and independent oscillations to both Wilson lines v* and v. Consequently 
the dipole amplitude s r (x — y) is expected to be small, thus enabling the non-linear term in 
( |2.11| ) to be dropped. The solution of the resulting linear equation can be shown to be M 



s T (x 



exp 



a K N r . 



dy\n(Q 2 s (y)x 



IT Jt 

If Q 2 s (y) = exp[ca s N c (y - r )/n}Q 2 (r ) pf, then 



(2-13) 



St X 



exp 



c fa s N c y , 
2{—) (r " ro) 



(2.14) 



provided r is chosen so that Qg(r )(x — y) 2 = 1. This condition implies (x — y) 2 = 
Q~ 2 (r ) ^> Qj 2 (t), if r 3> r . The solution (|2.14|) then confirms the expectation that s T is 
small in that region. 

To solve any of these evolution equations we need an initial condition at some r = r . 
For a large nucleus, and a r where the source is dominated by the valence quarks, the 
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initial condition is provided by the McLerran-Venugopalan model |Tj|, in which a Gaussian 
distribution is assumed for W T0 (a). For large A and small a s (Q 2 ), the Gaussian distribution 
can be shown to be a good approximation |17|. Saturation is now provided by the valence 



quarks alone, so A has to be very large and the resulting Q 2 S is relatively small. The detail 
of confinement which affects the small k 2 region then becomes relatively important [|TJ . 

With the Gaussian distribution, the perturbative gluon distribution is modified by a 
gluon dipole factor M 

^"W = j^^(^(x-,v)U(x-,0)), (2.15) 



where W is equal to in (|2.3| ) with r replaced by x , and r is the conjugate variable to the 



transverse momentum k. This factor gives rise to the non-linear effect that is responsible 
for saturation, with a saturation momentum Q s given by (|2.2|) . 

For r > To, the W T (a) determined by ( |2.4| ) no longer has a Gaussian distribution. Nev- 
ertheless, in a mean field approximation, the approximate solution is still Gaussian. The 
gluon distribution for small k 2 is then given by 

i(x r G)„ c2 Nt-l._K ln QM t ( k s«Q 2 ( T )) (2. 16 ) 



d 2 k * N c a s k 2 
for some constant c^- 



III. WILSON LINE DISTRIBUTION 

We shall show later in this section that the distribution W T (a) of YM potential a leads 
to a distribution W T (v ) of the Wilson lines. 

The notation might suggest that W T {y) describes only the distribution of anti-quark 
Wilson line v, but actually it provides a distribution for the Wilson lines of other partons 
as well. Since v is unitary, v* = t> _1 , the variable v' is a rational function of Vij, so W T {v) 
does provide the distribution for quark Wilson lines v\ 

From the group-theoretical relation 

t a Vl = vH c v = V cb t b) (3.1) 

or equivalently, the relation 

Vl = 2tv(vH c vt a ) = V ca , (3.2) 

gluon Wilson lines can be expressed as a quark anti-quark pair of Wilson lines, so we can 
also compute the distribution of gluon Wilson lines from W T (v). 

Physical observables are often expressible in terms of Wilson lines, so it is clearly desirable 
to know their distributions directly. Moreover, the Wilson line v is unitary, which allows 
the theory of representation of the unitary group to be used for computations. For example, 
an orthonormal complete set of polynomials, given by the irreducible representations of the 
unitary group, exists on the group manifold. Therefore a harmonic analysis of W T (v ) and 
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other physical quantities can be carried out to allow their integrals to be computed. In 
contrast, W T (a) is a function on a non-compact linear space of the Lie algebra, and the only 
functional that can be integrated in practice is the Gaussian. 

We shall show that W T (v) satisfies the same JIMWLK equation Q2.4D , but with H re- 
placed by H. The latter is obtained from H simply by replacing the a derivative by a 
differential operator in v. Therefore, H is still hermitean and positive semi-definite. 

Since the group manifold is compact, the spectrum and eigenf unctions of H are more 
manageable. For example, W T (v) = 1 is a normalized eigenfunction of H with zero eigen- 
value, whereas W T (a) = 1 has a divergent integral in the a space. 

Before showing how W T (v) is obtained from W T (a), let us first review some basic facts 
about integrations and orthonormal relations on a compact Lie group. 



A. Inner products on the color group 

For the sake of definiteness we shall assume the color group to be U(n), though a similar 
analysis can be carried out for SU(n). The n 2 -dimensional group U(n) will be parameterized 
by the n 2 matrix elements Vij (1 < i,j < n) in the defining representation. Unitarity equates 
to v~ x , so v*j is to be regarded as a dependent variable, given as a rational function of vm 
via this relation. 

The invariant volume element on U(n) will be denoted by e^M- This Haar measure is 
a left and right invariant n 2 -form, satisfying 

d H [v] = d H [v v] = d H [vv ] (3.3) 

for any constant vq G U(n). It is positive, and it shall be normalized to / dn[v] = 1. 

The Haar measure dn[v] is proportional but not equal to the product measure d[v], 
obtained by taking the exterior product of the n 2 1-forms dvij. They differ by a Jacobian 
J(v), so 

d H [v] = J(v)d[v]. (3.4) 

To get an idea how this comes about, consider a change dv in the vicinity of a group 
element v G U(n). Then v~ x dv constitutes a change around the identity, so it can be 
parametrized in the form —it a (dr] a ). The volume element at the identity is proportional to 
the exterior product of the n 2 drja's, or equivalently, the exterior product of 2iti(v~ 1 dv). 
If we want the volume element to be left and right invariant, as in (|3.3|) , this expression 
should be taken as the volume element at any point v G U(n). The presence of v~ x in this 



expression is the origin of the Jacobian J(y) in (3.4). When n is odd, dn[v] has a very simple 
analytical form, given by ( |A1[ ) in Appendix A. 

Let f(v) and g(v) be two functions on the group manifold. Their inner product is defined 
to be 

(f(v)\g(v)}= [ d H [v]f*(v)g(v). (3.5) 



It will be shown later that the operator D = Vij(d/dvij) is hermitean with respect to this in- 
ner product. Since Vvui = i>jw, the eigenf unctions of V are monomials of v , whose eigenvalues 
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are the degrees of the monomials. The inner product of two eigenfunctions of a hermitean 
operator is zero if their eigenvalues are different. Hence if M^iy) denotes a monomial of 
degree k and M[(v) a monomial of degree / 7^ k, then 

(M k (v)\Mftv)) = (k^l). (3.6) 

This result can be stated in another way. An integral / djj[v]B{v, v*) is non-zero only when 
the number of u's in B is equal to its number of v^s. We shall refer to this later as the 
matching rule. It is one of the main tools for our later calculations. 

To compute inner products when k = I, we resort to the theory of representation of the 
U(n) group, which asserts that if Dq D (v) is the matrix element of an irreducible represen- 
tation A, then 

J d H [v}Df B (v)D x CD (v) = J—S^SAcSsn, (3.7) 

where N[X] is the dimension of the irreducible representation A. 

For example, the defining representation v is irreducible and has dimension N c , hence 

(vij\vki) = jj-SikSji- ( 3 - 8 ) 

The orthonormal relation Q3.7|) can be used to compute inner products of any two mono- 
mials in the following way. First, apply Young's idempotent operators of k boxes to decom- 
pose Mk(v) into a linear combination of irreducible representations D x (v). Similarly M({v) 
is decomposed into a linear combination of D x (v). Then ( |3.7| ) can be used to calculate 
(Mk(v)\Mi(v)). Clearly it can also be used to calculate integrals / d H [v\B{y,v^). 



B. Measures and functionals on A and U 

To apply these properties to the physical problem on hand, we need to generalize them 
to the case when v depends on the transverse position x. We shall denote the color group 
containing v(x) as £/(n) x , and U = fix U(n) x . 

The product over x is to be interpreted in the following way. Cover the transverse x-plane 
by a square lattice with a lattice constant a. The product x is to be taken over all lattice 
points within the Lorentz-contracted nucleus of transverse radius Ra- The same convention 
will be applied to sums over x. 

The measure on U is defined to be 

V H [v]=l[d H [v{ic)), (3.9) 

X 

where djy[v(x)] is the Haar measure on Z7(n) x . Using ( |3.4| ) and denoting Ilx^(' y ( x )) by 
J(v), we get 

T>h[v] = J(v)V[v], 
V[v]=Y[d[v(jL)). (3.10) 
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The measure T>[a] on the Lie algebra A of YM potentials can be defined in the following 
way. Divide the y axis into intervals of size e. Since we are interested in a^J(x) only for 
y < t, the appropriate measure is 

V[a] =V r [a]V<[a], (3.11) 

where 

A- [a] =Y[da^(x), 

a,x 

D<[a] =c 3 I] II<K( X )> ( 3 - 12 ) 

y<r-e a,X 

with a normalization constant C3 to be chosen later. In the same vein, the Wilson line v (x) 
of (|2.3|) can be factorized into 



v(x) = t><(x)t> T (x), (3.13) 

with 

M X ) = ^ eX P y dy «y( X )^a) , 

u T (x) = exp (-i#<(x)t a e) . (3.14) 

A change of a"(x) of amount cfa;"(x) causes a change of f T (x) of amount efov(x) = 
— igev T ('x)t a da^('x) , and hence a change in i>(x) by the amount (if (x) = —igev('x)tada^('x). 
Or, f _1 <if(x) = —igeb a da,®. In light of the remark below ( |3.4j ), we can now choose the 
constant c 3 in ( |3.12| ) so that 

V H [v] = V T [a}. (3.15) 

We are now ready to discuss how W T (v ) can be obtained from W T (a). 

A functional of a y (x) for y < r can be folded into a functional of t>(x) using the formula 

F(v) = J F(a)5(v-u) jrtf>[a\ 



where 



and 



= jF(a)S{v- u)P[u]P<[a], (3.16) 
w(x) = Pexp (~i9 f dV , (3.17) 



<5(t; - u) = J] II 5 K x k " M ( x )^ > ( 3 - 18 ) 

x ij=l 

so that / f(u)6(v — u)V[u] = f(v) for any functional f(u). The second equality of (|3.16|) 
comes from flODD , ( pfi2|) , and ( gjg ). 
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It follows from ( |3.16|) that 

F{v)V H [v] = I V H [v]V[a]5{v - u)F[a]/J{u) 



F(a)V[a], (3.19) 

so W T (v ) is normalized if W T (a) is. 

It also follows from ( |3.16| ) by an integration by parts that the transform of SF(a)/Sa ( ^('x) 
is V a (yL)F(v), where 



£> a (x) = igvij(yi)(t a ) jk - — — = igti 
and the functional derivative in v is defined so that 



5v T (x) 



(3.20) 



5Uy(x) 



5piM(x-y). (3.21) 



The inner product between two functionals of v is defined to be 

(7(v)\g(v)) = jT(v)9(v)V H [v}. (3.22) 

It can be shown (see Appendix A) that £> a (x) is anti-hermitean with respect to this inner 
product. In particular, 



£>o(x) = ^•(x)[5/^(x)]/J2AT < 



is anti-hermitean, so the operator T> defined below eq. ( |3.5| ) is hermitean, as previously 
claimed. 

Instead of £> a (x) in ( ^20]) , it is more convenient to deal with the matrix operators f (x) 
and ^'(x), whose (mn) matrix elements are defined to be 

_ 2 — 5 

T-^mnip^) = ~ (ta)mnD a (x) fj n (x) — - -, 

ig dv im (x) 

c »>»< x > ^»<< x >fcir < 3 - 23) 

It can be checked that these two are related by 

t;(x)X)(x)t;t(x) = D / (x), (3.24) 

where the subscript x in v £(x) indicates that this v ^(x) should not be differentiated by the 
5/5v (x) in T>(x). In other words, in component forms, ( f3.24| ) reads fjmf^^mn = 
When operated on v pq and = v~ q x , they yield 

V mn {x)v pq (y) = +<5(x - y)5 mq v pn (y), 
V mn (x.)vl Q {y) = -<J(x - y)<V4 g (y), 
^mn( x )^(y) = +5(x - y)5 rv v mq (y), 

^Wi(y) = -*(x - y)<W<(y)- (3.25) 
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These formulas give rise to the following formulas which are very useful in practical calcu- 
lations. When D or T>' operates on a v or v* in a trace, we have 

D(x)tr [Av(y)} = +<5(x - y) [Av{y)\ , 
X>(x)tr M(y)A] = -5(x - y) [v\y)A 
V'(x)tr[v{y)A) = +5(x - y) [v(y)A] 



D'(x)tr Uv f (y)l = -5{x - y) \Av\y) 



(3.26) 



In other words, when the trace is written in a certain order, the operators -D(x) and -D'(x) 
simply remove the trace, and append to it the factor ±<5(x — y). When T> or T>' operates on 
a v or v ' in the same trace, we get 

tr [AV(x.)Bv(y)] = +5(x - y)tr [A] tr [Bv(y)) , 



tr \AV(^)Bv ] {y)\ = -5(x - y)tr [B] tr \v \y)A 
tr \AV'(x)Bv(y)} = +5(x - y)tr [B] tr [v(y)A] 
tr Up'(x) J Bi; t (y)l = -<J(x - y)tr [A] tr [Sv^y) 



(3.27) 



In other words, the single trace is broken up into a product of two traces. The matrices 
A, B in these equations are constant matrices. 



C. W T (v) and its JIMWLK equation 



Using (|3.16|) , the distribution function W T (a) can be folded into the distribution function 



W T (v) of Wilson lines. Since <J/<fa£(x) is transformed into D a (x) of (^20|) , the JIMWLK 
equation ( |2.4| ) for W T (a) is now changed into a JIMWLK equation for W T (v): 

d ^ = -HW A v), (3.28) 

where 

H= - I d 2 zd 2 xd 2 y^(xy|z)0(xy|z), 

7T J 

0(xy|z) = D a (x) (yt( x ) _ (z) )^ (y(y) _ v(z))cfe D b (y). 

(3.29) 

Like if, if is also hermitean and positive semi-definite. 

Using ( |3.23| ) and ( |3.2|) , the operator 0(xy|z) can be written in a form more convenient 



for practical calculations. Remember for this purpose the remark after ( |2.7|) , that the factor 
[Vt(x) — V"t(z)] ac in ( |2~5| ) can be written to the left of the differential operator <5/<5a"(x). 



In terms of (|3.29|) , this means to the left of 7J> a (x). In what follows we shall use 0(x.y\z) of 
( |3.29| ) in this form. 

Using (|3.23|) , (|3.24|) and ( |3.2|) , and this remark above, we have 
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V cb (z)V b (y) = igti (V(z)M(z)X?(y) 
V a {y)V\z) = igti (v\z)t c v(z)V(x)) , 



V cb (y)V b (y) = igtr (v\y)t e v(y)V(yj) = igti (t c V'(y)) , 
P a (x)y f (x) = igti (v\x)t c v(x)V(x)) = igtr{t c V'{x)) . 

From these relations, and the identity ti(t c A)tr(t c B) = ^tr(AB), we get 



0(xy|z 
O 



o xy + o xz + o zy + o 



zy 



■HI 



--ghr[D'(x)V'(y)}, 
+~<? 2 tr [&{x)v(z)V{y)vfa) 
O yz = +^ 2 tr U(z)2?(x)4(z)P'(y) 



g 2 tr v(z)V(x)vl(z)v(z)V(y)vUz 



-^ 2 tr[P(x)P(y)]. 



(3.30) 



(3.31) 



Assuming W T {v) to be normalized, JVh[v]W t (v) = 1, the average of any functional 
B(v, v') is equal to 



(B) T = J V H [v]B(v,v^)W T (v). 



(3.32) 



If B and W T are both monomial functionals of v and v', this functional integral factorizes 
into a product of integrals on the group U(n), each of which can be computed using ( |3.6|) 
and ( |3.7| ). In particular, the functional integral is non-zero only when the number of v's and 
uT's in BW T are the same for every transverse position x. This is the functional form of the 
matching rule previously considered. 



IV. MUELLER'S DERIVATION OF THE JIMWLK EQUATION 

In a recent paper |J, Mueller proposed a simple derivation of the JIMWLK equation in 
the following way. He showed that if W T (a) is equal to the dipole functional v^(s)ijV (t) M , 
then the Feynman diagrams for dW T (a)/dr can be written in the form —HW r (a), with H 
given by ( |2.5|) . He then stated that the same is true if W T (a) is equal to any multipole 
functional (si)^^ X s 2)i 2 j2 ' ' ' v i^i)k 1 i 1 v{t2)k 2 i2 ' ' ■> an d hence it is likely that the JIMWLK 
equation for W T (a) is also valid. 

In this section we shall supply the missing steps of this proof. This consists of filling 
in the detailed arguments for the multiple functionals, and then showing that they lead to 
the JIMWLK equation for W T (v). We do not know whether the JIMWLK equation for 
W T (a) follows or not. However, for all our applications, a JIMWLK equation for W T (v) is 
sufficient, so it really does not matter whether the equation for W T (a) can be proven this 
way or not. 
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Instead of starting from the Feynman diagrams to derive the evolution equation, as is 
done in Ref. ||, we find it easier to do everything in reverse. That is, we start from the 
JIMWLK equation and show that they lead to the correct set of Feynman graphs. This 
inverse approach makes it more manageable to deal with the complicated multipole func- 
tionals. Actually, simplification already occurs at the dipole level: a necessary cancelation 
in the original derivation is avoided altogether in this way. 

The Wilson line v'(s) for a quark is drawn in Fig. 1 as a left-pointing arrow, and the 
Wilson line v(t) for an anti-quark is drawn as a right-pointing arrow. Time x~ is drawn 
to increase from right to left; multiplication of color matrices from left to right should be 
carried out against the arrow of the fermion lines. 

The short vertical bars in the middle of the lines (labelled IP) indicate the lightcone 
longitudinal position x~ = where interaction between the multipole and the pancake 
nucleus takes place. Since a(x) is concentrated around x~ = 0, we may regard v and to 
be located at the IP. 

The operation of D a (x) or 5/5a®(x) on the Wilson lines is given in ( |2.7|) . This operation 
can be represented graphically by putting a vertex to the left (the larger x~ side) of IP, both 
for the quark Wilson line and the anti-quark Wilson line v. The vertex for a quark is 
igt a , and the vertex for an anti-quark is —igt a . 



Using the remark following eq. ( |2.7| ), the operator 0(xy|z) in (2J3) can be written in the 
form 

=[^(x)-^(.)] - ^.[7(y)-V(.)U^ y 

= O xy + O xz + O zy + O zz . (4.1) 

Using ( ft.lD ) and the unitarity relation vv^ = v^v = 1, we see that the effect of 
[U(y)] c& (S/Sa^.(y)), operating on a Wilson line, is to place a vertex to the right (the small 
x~ side) of the IP on the Wilson line. Similarly, the effect of V^(x) (5/5a®(x)), operating 
on a Wilson line, is also to place a vertex to the right (the small x~ side) of the IP on the 
Wilson line. 

We are now ready to see what happens when 

H = J d 2 zd 2 xd 2 yK(xy\z) (O xy + O xz + O zy + O zz ) 

= H xy + H xz + H zy + H zz (4.2) 

operates on a multiple functional, i.e., a collection of p quark and q anti-quark Wilson lines. 

H xy puts a vertex y to the right of the IP on a Wilson line, and vertex x also to the right 
of the IP of the same or a different Wilson line. These two vertices are then connected by 
the 'gluon propagator' / d 2 zK(-xy\z)5 a b, where a, b are the color indices at the two vertices. 
The gluon propagator is shown in Fig. 1(a) with a dashed line. This operation is to be 
applied to every pair of Wilson lines, including the possibility of applying to the same line 
twice. 

Similarly, since Vt(z)V(z) = 1, the effect of H zz is to put a vertex y to the left of the IP 
of a Wilson line, and another vertex x to the left of the IP of the same or another Wilson 
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line. The two vertices are then linked by the same 'gluon propagator'. This is shown in 
Fig. 1(b). 

In both of these cases, the two vertices are both to the same side of the IP. This is not 
the case with the other two terms. 

H xz puts a vertex y with color b to the left of the IP, and a vertex x with color a to the 
right of the IP. These two vertices are linked by a 'gluon propagator' / d 2 z [V(z)] ab K(x.y\z). 
This is shown in Fig. 1(c). 

There is potentially another contribution to if xz when V^ c (x)6 / 5a a (x) operates on V c b(z). 
However, this term is proportional to <5(x — z)V^ c (x)V c d(z)(T a ) db , which is proportional to 
{Tajab = 0, so that term is actually absent. 

Similarly, H zy puts a vertex y of color b to the right of the IP, and a vertex x with 
color a to the left of the IP. These two vertices are linked by the 'gluon propagator' 
/ d 2 z V'(z) K(xy\z). This is shown in Fig. 1(d). 

L 1 ab 

Appropriate signs and Dirac ^-functions on the transverse coordinates must also be 
inserted. 

One might be bothered that the gluon propagators in the four diagrams appear to be 
different. Fortunately this is only superficial. To see why they are actually the same, first 
note that T a is imaginary so V is real. Hence the gluon Wilson line [V^(z)] a & in the propagator 
in Fig. 1(d) can be written as [V(z)]& a . If we compare this with that of Fig. 1(c), we see 
that these two are the same, both equal to [V(z)] 0l02 , where a\ is the color index before the 
interaction point IP, and a<i is the color index after. Since a a (x~,z) is concentrated near 
x~ = 0, we may write, in both cases, the gluon Wilson line to be 



[^( z )ka 2 = ^Pexp 



f x 2 

—ig / dx~a a (x~, z)T a 

Jx7 



(4.3) 



where (ai,xj~) is the interaction vertex to the right of IP, and (02,^2") is the one to the left 
of IP. For the first two diagrams, we may replace 5 a b by the same expression (|4.3|) , because 
in that case both Xi and x^ are to the same side of IP, hence a a (x~, z) = throughout the 
integration interval, so [V(z)] & = 6 a b- 

These four types of Feynman graphs are precisely those needed for the development 
dW T (a)/dr @. Therefore the JIMWLK equation is satisfied whenever W T (a) is given by a 
multipole function, namely, a monomial functional of v and v\ Since polynomial functionals 
on the group manifold form an orthonormal complete set, the JIMWLK equation ( 3.28|) 



must be true in general. This completes Mueller's proof for W T {v). However, since W T {a) 
depends on many more variables a a (x~,x) than f(x), it does not necessarily follow from 
this argument that (|2.4j) for W T (a) is true. 



V. INFRARED DIVERGENCE AND THE GENERALIZED BK EQUATION 

The kernel K(xy\z) of the JIMWLK equation, given in ( |2.6|) , goes like 1/z 2 for large 
|z|. This causes a log divergence in HW T (a) of fl2.5|) , and H W T (v) of ( |3.29| ). On the other 
hand, the dipole amplitude ( |6.1|) satisfies the Balitsky equation fl2.10|) , whose kernel goes 
like l/(z 2 ) 2 for large |z|, so infrared divergence is absent in that case. In this section, we 
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shall use the Feynman diagrams derived in the last section to show that the multipole traces 
defined below are equally free of infrared divergence. 
Multipole traces are defined by 

m(siti • ■ ■ s fc t A ) = — tr |y(si)u(ti) ■ ■ ■ v^(s k )v(t k ) . (5.1) 

We will call the functional average of a multiple trace, (mjt(siti • • -s k t k )), a multiple ampli- 
tude. When k = 1, this reduces to the dipole amplitude Q2.8|) . 

The diagram for a multiple trace is shown in Fig. 2. As far as color-matrix multiplication 
is concerned, these 2k lines should be considered to be joined together at their ends to form 
a single big loop with the arrows all pointing the same way. The joints are indicated by 
dotted lines in the figure. Note that if the amplitude did not have the form displayed in 
flS.ip, with and v occuring alternately inside a trace, such a big loop cannot be formed 
and the argument below will not be valid. 

Consider the two gluon propagators shown in Fig. 2. The color structures are identical, 
but the two terms differ by a sign because the vertex B is on a quark line in one diagram 
and on an anti-quark line in another diagram. Their gluon propagators may therefore be 
combined into K(xx.\z) — K(x.y\z). For large |z|, this is proportional to z-(x — y)/(z 2 ) 2 = 
0(l/|z| 3 ), hence the infrared divergence is absent. As a matter of fact, we can also combine 
diagrams with vertex A similarly shifted. The four K thus combined actually dies down like 
l/(z 2 ) 2 for large |z|. 

It is clear that all the diagrams for the multipole function can be paired up in a way 
similar to Fig. 2, thus eliminating all infrared divergences. Moreover, since the combination 
is obtained by combining two diagrams with vertex B in different positions, but the same 
vertex A, finiteness persists for products of multipole traces. 

In the rest of this section, we shall derive a generalization of the Balitsky-Kovchegov 
(BK) equation to multipole amplitudes. We will see that the kernel of the equation actually 
behaves like l/(z 2 ) 2 for large |z|, just like the kernel for the Balitsky equation. 

Since H is hermitean, the JIMWLK equation satisfied by W T {y) is also satisfied by 
the multipole trace flOP . We may now use (Q), (|3T2p , (|3~29|) and ( $M ) to simplify 
Hm{siti ■ ■ ■ Sfct/j). The resulting equation is 

<9m(siti ■ • -Sfctfe) 



dr 



H xy + H zz + H xz + H yz ] m(sita ■ ■ • s k t k ), (5.2) 



where 



1 ^ 

H xy m{ Sl ■ ■ ■ t fc ) = +^9 2 N C E / d 2 zl ij {z)m^{s l ■ ■ ■ t k ) 

1 ^ 

H zz m{ Sl ■ ■ ■ t fc ) = +7>9 2 N C E / d 2 z^(z)mg( Sl • • -t k ) 



1 ^ 

H xz m{ Sl ■ ■ ■ tfc) = - -g 2 N c E / d 2 zl lj {z)^{s 1 ■ ■ ■ t k ) 



1 ^ 

H zy m{ Sl ■ • • tfc) = - -g 2 N c E / d 2 zl lj {z)m^{s 1 ■ ■ ■ t k ). (5.3) 



16 



The kernel in these equations is 



Iij(z) = K(siSj\x) - K(s i t j \z) - K(sjti\z) + K(titj\z), (5.4) 

which is 0(l/(z 2 ) 2 ) for large |z|. It is also symmetric in % and j: hj{z) = Iji(z). The 
meaning of the multipole traces in (|5.2| ) will now be explained. 

The argument (s! • • • t fc ) of the multiple trace in ( |5.1| ) is circular, because the trace is. We 
can consider the argument s x to be behind the argument t^, forming a circle. For example, 

m^tx ■ • -Sfetfc) = m(tis 2 • • -t fc Si) = -^rtr ^(ti)^^) • • •u(t fe )u t (si)] . (5.5) 

In these formulas, v is always associated with a tj and v* is always associated with a s$. 

The quantity m^ a (si • ■ • t*.) appearing in the if xy term in ( ^3|) is defined as follows. Put 
a vertical bar after Sj, and another vertical bar after Sj in the circular argument (si • • • t^). 
This pair of vertical bars separates the circular argument into two circular arguments. 
m^(si ■ ■ ■ tfc) is defined to be the product of two multiple traces with these two circular 
argument. For example, 

m"3(s 1 t 1 s 2 t2S3t 3 ) = m(si|tiS2t 2 S3|t 3 ) = m(t 3 si)m(t 1 s 2 t2S3). (5.6) 

The superscript a stands for 'after'. The superscript b to be found in the other three 



terms of (|5.3|) stands for 'before'. Each superscript together with its corresponding subscript 
pair up to tell us where the vertical bar is put. The pair (ai) means to put a vertical bar 
after Sj. The pair (bj) tells us to put a vertical bar before Sj. With this understanding the 
quantity ?n^(si • • -t^) can be defined similarly. For example, 

m5 3 (sitis 2 t 2 s 3 t 3 ) = m(|sitis 2 t 2 |s 3 t 3 ) = m(sitis 2 t 2 )ra(s 3 t 3 ). (5.7) 

If we apply these recipes literally to the last two equations in (|5.3| ) we will end up with 
something non-sensical, in that the numbers of v's and 's within each trace is not identical, 
thus neither is a multipole trace. What we should do in these two cases is to insert a v(z) into 
one trace and a v\z) into another so that both become multiple functionals. For example, 

m" 3 (sitis 2 t 2 s 3 t 3 ) = m(s 1 |t 1 s 2 t 2 |s 3 t 3 ) = m(t 1 s 2 t 2 z)m(zs 3 t 3 s 1 ). (5.8) 

The z in the first factor is a v'(z), and that in the second factor is a n(z). 

This completes the explanation of the symbols in ( p.3[ ), except for one last remark. If 
the two vertical bars are side by side, then the corresponding mulitple trace m should be 
interpreted as tr(l)/iV c = 1. 

When we take expectation values on both sides of (|5.2| ), we get an equation for the 
multiple amplitudes. For large N c , the expectation of products of traces factorizes into 
products of expectations of a trace. In this form ( [5.2| ) and ( |5.3| ) generalize the Balitsky- 
Kovchegov (BK) equation to higher multiple amplitudes. We should interpret m in these 
equations as the average (m), and products of two m's as the product of the averages. 

In the special case of a dipole, k — 1, hence i — j — 1 in (|5.3|) . In that case, 
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1 1 


= m(si ti) 

\ -Ml - 1 - / 




= mdlsiti) 




= m c fl = m 




1 ( 

47T 3 (Si - 



m(siti), 
m(siti), 

i\ti) = m(s 1 z)m(zt 1 ) 



Sl -t ^ 2 



-. (5.9) 

Equation (|5.2|) coincides with the BK equation ( p.ll| ), as it should. 

A very interesting fact emerges from these generalized BK equations for multipole am- 
plitudes. If the multipole amplitudes for k = 1,2, ■■■,£— 1 are known, then the evolution 
equation determining the £th multipole is a linear inhomogeneous equation. Thus the only 
non-linear equation one has to solve is the original BK equation for the dipole amplitude. 
It is well known that the non-linearity of this equation leads to saturation. The linear- 
ity of the higher multipole amplitudes therefore means saturation of the dipole amplitude 
automatically drives saturation of all higher multipole amplitudes. 

If we require the solution of the JIMWLK equation ( |3.28j ) to be free of infrared singular- 
ities, then presumably it will be made out of the multiple traces (|5.1| ) and their products. 
In that case the difficult functional differential JIMWLK equation can be replaced by the 
more manageable set of generalized BK equations (|5.2j) and (|5.3|). 



VI. SOLUTION OF THE JIMWLK EQUATION 

Suppose we decompose W T (v) into a linear combination of eigenfunctions <fi\(v) of H. 
The operator H is hermitean and positive semi-definite, hence the eigenvalues A are real and 
non- negative. From (|3.28|) , the r dependence of W T (v ) is thus given by a linear combination 



of exp(— Ar). In the asymptotic limit r — > oo, the lowest eigenvalue of H dominates. 

The lowest eigenvalue is A = 0, and its normalized eigenfunction is 4>o(v) = 1. This 
eigenfunction is normalized because / T>h[v] = I\ X J dn[v (x)] = 1. 

In contrast, 0o( a ) = 1 is also an eigenfunction of H with A = 0, but the integral of this 
eigenfunction is divergent because the linear space A is non-compact. 

In the next subsection, we will discuss what happens in the asymtotic limit r — >• oo when 
A = dominates. In the subsequent subsection, we will look at W T (v) and its applications 
for a smaller r. 



A. Asymptotic limit 



When r = oo, only 



contributes, so we can set W r 



1. Averages 



are then given by the integral (B(v } v*)} = J T>h[v }B(v, v '). From the matching rule, this 
integral is non-zero only when the number of w's in B exactly matches the number of tr's 
at every transverse position x. 

For example, the dipole amplitude (|2.8j) is 



«r(x,y) = — (tr [v t (x)«(y)J) = 5 x , y 



(6.1) 
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because v^v x = 1 and tr(l) = N c . This is consistent with (|2.9|) , and also ( p.!4|) at r = oo. 
In other words, whenever the dipole has a finite size, the dense gluon will have such a strong 
absorption that the dipole amplitude always become zero. 

A similar statement can be made about multipole amplitudes. 

Let us now compute the gluon spectrum at r = oo by using the formula B 



d(x F G) 



d 2 k 4vr 3 
where k = (fc + ,k), and 



1 



4tt 3 N c 



(Tr F + \k)F + \-k) ) 



F+\k) = / rf 3 xexp(ik-x)F 



X 



(6.2) 



(6.3) 



is the color electric field in the lightcone gauge. In the approximation F +l (x) ~ 
(i / g)5(x~)V (x)d l V^ (x) which is supposed to be valid for a Lorentz-contracted pancake nu- 
cleon, the formula becomes 

d(x F G) 



d 2 k 
C(x,y) 



J d 2 xd 2 yexp[ik-(x- y)](C(x,y)), 



Tr 



V(x)V?(x)V(y)V?(y) 



(6.4) 



where the subscript i represents a differentiation, namely, V} = diV^ . 

To compute the average (C(x, y)) = / T>h\p\B{x., y), we need to use 
gluon Wilson lines V and to the quark and anti-quark Wilson lines v' and v. The result 



to convert the 



is 



C(x,y) = 2C F tr \v(y)v}(y)v(x.)vj(x) + v t (y)v\y)vi(x)v\x) 



+ tr 
+ tr 



tr 
tr 



u(x)v}(x 
Vi(x.)vUx 



The functional integral (C) = / T>h[v\C is computed in Appendix B. The result is 



<C(x,y)) = -2 



N 2 



-5- 



(6.5) 
(6.6) 



Changing the integral in 



d(x F G) 
d 2 k 



into a finite sum, and letting r = x — y, we get 

= W^? aexp(2k ' r) ^^^ 

2vr J R 2 4 (iV 2 - 1) 



16vr 4 a s iV c ' ^ 

The unintegrated spectrum d(x F G)/dk 2 is therefore absolutely flat, up to the saturation 
momentum Q 2 which is in this case infinite. That is not unexpected at r = oo. 
The integrated cross section is then 



x F G = Q 



2 ttR\(N 2 -1) 



(6.8) 



The spectrum in ( |6.7| ) differs from the mean field prediction ( |2.16| ), in that ( |6.7| ) is flat 
and ( p,16| ) has a logarithmic dependence on k 2 . The integrated density x F G is however quite 
similar to the estimate given in (j2.2|); the only difference is a factor c = 8ir 3 . 
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B. Below the asymptotic limit 



It is much more difficult to solve the JIMWLK equation for finite r, because we know 
nothing about the other eigenfunctions and eigenvalues of H. In this section, we will discuss 
an approximate correction to W T (v) below the asymptotic limit, in a region where the 
Wilson lines are far apart. 

To avoid the infrared divergence, the distribution functional W T (v) will be assumed 
to depend on v and v' only through the multipole traces ( |5.1| ), or products of them. From 
results of the last subsection, and discussions in Sec. 2, we know that for large r, the average 
of multipole traces (i.e., multipole amplitudes) are small if the Wilson lines in the multipoles 
are far apart. In fact, the higher the order of the multipole is, the smaller the amplitude is 
expected to be. Therefore it is reasonable to include only quadratic dependences of v and 
v* in a first correction to the asymptotic limit, at least in the region when the Wilson lines 
are far apart. We will therefore assume 



s.t 



s)v(t) 



6r(t,s). 



(6.9) 



The contribution tr[v t(s)u (s)]6 T (s, s) 



N c b T (s,s) may be absorbed into the 



independent term. We may therefore assume 6 T (s, s) = without any loss of generality. 
In that case, using the matching rule and unitarity of v to do the functional integral, we see 
that W T {v) is still normalized: 



J V H [v\ W T (v) 



1. 



(6.10) 



To compute the dipole amplitude (|2.8|) , we need the following integration formula, which 
can be obtained from the matching rule, eq. (|3.8|), and the unitarity of v: 



V H [v] tr 



vHs)v(t) 



tr 



f (x)t;(y)l = iV c 2 5 Sit <5 x , y + 5 s , y 5 t ,x - S t 



s,y,t,x- 



(6.11) 



The last Kronecker 5 is by definition non-zero only when the four arguments in its subscript 
are all equal. 

We may now compute the dipole amplitude from (|2.8| ) to be 



M x >y) 



Nr. 



V H [v] W T (v)ti u t (x)v(y) 



^x, y + — 6 T (x,y) 



(6.12) 



Since 6 r (x, x) = 0, we get s T (x, x) = 1, as it should be. If we stay away from x = y, 
then 6 r (x, y)/N c is just the dipole amplitude s T (x, y). As such it should satisfy the BK 
equation. On the other hand, we should be able to get the equation of 6 T (x, y) directly 
from the JIMWLK equation ( g2§) by using ([OP , ( ggg ), and ( ggg ). The left hand side is 
proportional to tr v^(s)v(t) , but the right hand side has two terms, one is proportional to 
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tr 



vHs)v(t) 



\v{t) 



tr 



vHs)v(z) 



If we drop this 



and the other it proportional to tr 
quartic term because it involves a higher order multipole function which is expected to be 
small, then 6 T (x, y)/N c = s T (x, y) simply satisfies the BK equation ( [2.11| ) with the quadratic 
term of s T dropped. This is justified when the dipole amplitude is small, which is the case 
when the two Wilson lines are far apart, as expected. The solution is given by ( |2.13| ). 

We may use ( |6.9|) to calculate higher-order amplitude. The result is a sum of b T /N c = s T , 
one for each dipole pairs inside the multipole. 



VII. CONCLUSION 



Density of soft gluons is determined by the distribution W T (a) of the classical Yang- 
Mills potential a(x). Interaction of fast partons with such a background is given by their 
Wilson lines. In this paper, we introduced the distribution W T (v ) of Wilson lines. It can be 
obtained from W T (a), and it also satisfies the JIMWLK equation. We completed Mueller's 
derivation of the JIMWLK equation, though for W T {v) and not for W T (a). We derived a 
generalized BK equation for multipole amplitudes. We also used the normalizable property 
of W T (v ) to compute the properties of physical observables at x F = 0. We obtained in this 
way an unintegrated gluon spectrum d(x F G)/dk 2 = nR\(N 2 — l)/8ii 3 a s N c , independent 
of its transverse momentum k 2 of the gluon. Correction to this asymptotic behavior was 
briefly discussed. 
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APPENDIX A: MATHEMATICAL DETAILS 



The Haar measure of U(n) for odd n is given by 

dji[v] = C4 tr v~ x dv A v ~ dv A • • • A v _1 dv , (Al) 

where C4 is a normalization constant determined by / du[v] = 1, and the argument of tr 
consists of the exterior product of n 2 factors of v~ x dv. From the cyclic nature of the trace, 
and the antisymmetry character of the exterior product, it can be seen that dn[v] defined 
this way is when n is even. Hence we must confine ourselves to an odd n of we want to 
use the expression (|A1|) . 

It is easy to check that ( ^.3[ ) is satisfied for ( |A1| ). 

We will show in two ways that T> a (x) defined in ( |3.20| ) is anti-hermitean. First, define 
the inner product of two functionals a by 
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(f\g) Q = / V[a\V[a'\r{a)g{a')8 (u - v!) /J(u), (A2) 



where u is given by ( |3.17| ) and v! is given similarly in terms of a'. From ( p . 1 0| ) and fl3.15| ), 
we see that 

S(u - u')/J{u) = J] 5 (a T (x) - o4(x)) . (A3) 

Using integration by parts and assuming the resulting surface term to be zero, it is easy to 
see that 8/5a.®(r, y) is anti-hermitean with respect to this inner product. 

If f(v) and g(v) are obtained respectively from f(a) and g(a) by ( |3.16| ), then the inner 
product between f{v) and g(v) defined by 



(f\g) v = J V H [v]f(v)g(v) (A4) 

is equal to the inner product (f\g) a defined in ( |A2|) . This is so because 

(7\g)v = i V H [v]V[a\V[a']r(a)5(v - u)g(a)5(v - u')/J{u)J{u') 

' V[a]V[a']5(u - u')f{a)g{a')/ J{u) = (f\g) a . (A5) 



We know that the transform of 5f(a)/5a®(y) is t> a (y). Since (fg) v = (f\g)a and since 
5/5a®(y) is anti-hermitean with respect to (f\g) a , h follows that T> a (y) must also be anti- 
hermitean with respect to (f\g). 

The second proof of the anti-hermiticity of ^(y) makes use of the explicit formula ( |A1[) . 
It proceeds as follows. 

Using the explicit formula (|3.20| ) for V a (y), we must show that T> a (y) = 



v ij(y)(t a ) jkfi / dvik(y) is hermitean with respect to the inner product (|A4j) . 



Let us first show that T> a (y) is imaginary. Since v'—v 1 , it follows that 

v is v rk - — . (A6) 



5v* k Su 
Using also the fact that t a is hermitean, then 



Now we use integration by parts to show that T> a is antisymmetric with respect to the inner 
product (f\g) v . Remembering ( |3.10| ), integration by parts changes V a into 



J(v) 5v ik 

This would indeed be equal to —V a if 

8 



1 s 

(vt a ) ik J(v). (A8) 



~ r [(vt a ) ik J(v)\ = 0. (A9) 
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Since 5 (d[v\) /Svik = d [8v/8vn-\ = 0, ( |A~9|) is true if 



S [{yt a ) ik V H [v} 



8v 



ik 



which in turn is true if 



ntr(t a )d H [v] + (vta)*-^- = 0. 



dvi 



(A10) 



(All) 



ik 



Using (|A1|) , we get 



ik ' 



dd H [v] 
dv 



—n 2 tr 



ik 



t a v 1 dv A ■ • ■ A v 1 dv 



(A12) 



namely, it is equal to — n 2 times a t a inserted into the measure in (|AT|) . If a = 0, then to is 
proportional to the identity matrix, so indeed ( All ) is true. For a > 0, tr(£ ) = 0. Since 
v~ x dv is a linear combination of the U(n) generators tf,, we conclude that (|A11|) is valid if 



0. 



(A13) 



This is indeed the case because tr(t a ) = and s = €b 1 b 2 --b n2 tatb 1 • • • h 2 * s proportional to the 
identity matrix. This last statement can be proven as follows. 

Let Vq G U{n). Then v q t a v o = (Vo)a6^6, where V is the adjoint representation of v . 
Hence v^svo = det(Vo)s. Since the adjoint generator T a is has no diagonal matrix elements, 
it is traceless, so det(Vo) = 1. Therefore s commutes with every element v of U(n), so by 
Schur's lemma, it must be proportional to the identity matrix. 

We have thus shown that T> a (y) is antisymmetric and imaginary. Hence it is hermitean. 



APPENDIX B: SATURATION CALCULATION 

We want to compute (C(xy)) = /X>#[V]C(x, y) for the function C(xy) given in QS.5|). 
This function contains four terms. We will label them consecutively as Ci, C2, C3, and C4. 
We interpret the derivative div(x) = fj(x) on the lattice to be 

u<(x) = — [«(x + a<) -v(x-a<)], (Bl) 

where aj is a lattice vector along the ith direction, v j(x) is defined similarly. The following 
computations make use of (|3.6|) , Cf = N c /2 for U(n), the matching rule discussed at the 
end of Sec. 3, and unitarity of v. There are also extra factors of 2 obtained by summing the 
two aj directions. In this way we get 

a,)) 
a,))]) 



(Ci) 



4^2 



(tr 



>(y) (^(y + a-i) -v\y 



vix) [v'hc + a^ — ir (x 
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(C 2 > 



(C 3 > 



[Ox,y-a< + <W+aJ ^ -5— 3 



4a 2 
AT 2 

c r 

,2 ° x -y 



a 
2C F 

4^ 



(tr 



(u(y + a*) -w(y-a i ))w t (y) 
(u(x + a*) - u(x - a*)) u f (x) 



(Ci), 
1 



4a 2 



tr 



a 



tr 



(ft) = J^r 



tr 



(v(y + ai) - u(y - a*)) u f (y) 
u(x) (V(x + 3i) - v f (x - a,)) 

u(y) (^ f (y + a i) - vt (y - a *)) 

(u(x + Si) 



ViX — 8Li)V' (x 



(C(x, y)) = (d + C 2 + C 3 + C 4 > = -2 



JV?-1 
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FIGURES 
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FIG. 1. Diagrams representing the result of H = H xy + H zz + H xz + H zy , or 
H = H xy + H zz + H xz + H zy , operating on a multipole functional with 2 anti-quark Wilson 
lines (left-pointing arrows) and three quark Wilson lines v (right-pointing arrows). The dashed 
lines are the 'gluon propagators'. See the text in Sec. 4. for further explanation. 
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FIG. 2. A multiple trace with k = 3. The two diagrams have identical color structure, so their 
two 'gluon propagators' can be combined into an infrared finite expression. 
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